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ABSTRACT 

Na I D lines in the spectrum of the young binary KH 15D have been analyzed 
in detail. We find an excess absorption component that may be attributed to 
foreground interstellar absorption, and to gas possibly associated with the solids 
in the circumbinary disk. The derived column density is log A^ai = 12.5 cm" 2 , 
centered on a radial velocity that is consistent with the systemic velocity. Sub- 
tracting the likely contribution of the ISM leaves log A^Nai ~ 12.3 cm" 2 . There is 
no detectable change in the gas column density across the "knife edge" formed 
by the opaque grain disk, indicating that the gas and solids have very different 
scale heights, with the solids being highly settled. Our data support a picture of 
this circumbinary disk as being composed of a very thin particulate grain layer 
composed of millimeter-sized or larger objects that are settled within whatever 
remaining gas may be present. This phase of disk evolution has been hypothe- 
sized to exist as a prelude to the formation of planetesimals through gravitational 
fragmentation, and is expected to be short-lived if much gas were still present 
in such a disk. Our analysis also reveals the presence of excess Na I emission 
relative to the comparison spectrum at the radial velocity of the currently visible 
star that plausibly arises within the magnetosphere of this still-accreting young 
star. 
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Introduction 



KH 15D is a system with u nique and dramatic photometric behavior that was discov- 
ered by iKearns fe Herbstl ( 119981 ) during a variability study of the extremely young cluster 
NGC 2264. Basic properties of the star, including i ts distance (760 pc) an d classification as 
a weak-lined T Tauri star (WTTS) were reported by Hamilton et al.l (120011 ) . It is now known 
to be a young (~3 Myr) eccentric binary system embedded in a nearly edge- on circumbi 



nary disk that is tilted, warped and precessing with respect to the binary orbit (IWinn et al. 



2004 : IJohnson fc Winnll2004J : IChiang & Murrav-Clavil2004J : I Johnson et al.ll2004h . When first 
discovered in 1995, the disk was oriented relative to our line of sight in such a way that it 
covered most of the orbit of one star (now known as Star B) but only a small portion of the 
orbit of the other (Star A). During the past 14 years, precession of the disk has caused the 
sharp edge of the opaque screen to gradually cover all of the orbit of Star B and most of 
the orbit of Star A. As a result, we have not seen Star B (except by reflected light) since 
1995, and have seen direct light from Star A fo r progressive l y less of its 48.4 day orbital cycle 
([Hamilton et al.ll2005l ; iLeDuc et al.llin prep.1 ). IWinn et al.l (120061 ) have produced a model of 
the system that accurately accounts for most of the photometric and radial velocity obser- 
vations. They estimate the masses of Star A and Star B to be 0.6 and 0.7 and the radii 



to be 1.3 and 1.4 R^, respectively. The spectral class of Star A is K6/K7 (IHamilton et al. 
200ll : lAgol etUl2004h . 



The system's geometry and fortunate orientation have provided astronomers with a 
unique opportunity during the past decade and a half to stud y the structure of b oth the 
gas and the dust components of an evolved protoplanetary disk. iHerbst et al.l (120081 ) argued 
that the particulate grains making up the solid component of the disk, the "grain disk" , are 
probably of sand (millimeter) size or larger. It is possible that KH 15D possesses ex actly 
the sort of settled p articu late disk originally envisioned by ISafronov fc Zvjaginal (119691 ) and 
Goldreich fc Ward! (119731 ) as a prelude to the formation of planetesimals through gravita- 
tional instability (see also lYoudinl 120081 ) . The gas disk at this stage is predicted to have a 
much larger scale height than the particulate disk since it is partially supported by pressure. 
Because Star A of the KH 15D system changes elevations with respect to the edge of the 
grain disk, during the past decade we were presented with an opportunity to observe differ- 
ent path lengths through any accompanying gas disk for signs of absorption. We know that 
there is some gas still pres ent in the disk because the star has an activ e magnetosphere and 



drives a weak bipolar iet (IHamilton et al.ll2003l ; iTokunaga et al.l 12004c iDeming et al.l 12004 



Mundt et al.ll2010l : IHamilton et al.llin prep.l ). KH 15D has the potential to tell us about an 
important missing link in the planet formation process: the transition from protoplanetary 
disk to debris disk. 



-3 - 



In this paper, we constrain properties of the gas component of the KH 15D circumbinary 
disk by carefully analyzing the Na I D absorption lines visible in the spectrum of Star A while 
at various heights above the occulting disk. In §2] we describe the echelle spectra available for 
this analysis, their reduction, and the basic features of the observed Na I absorption profile. 
In §3.11 we discuss how Na I column densities are computed. In §3.21 we measure interstellar 
Na I absorption toward KH 15D and discuss how this affects our absorption measurements. 
In §4. II we discuss the implications for the scale height of the gas disk, and in §4. 21 we further 
discuss the interpretation of the data in terms of what it may reveal about the putative gas 
disk of KH 15D. Finally, in §5] we summarize the insight into the KH 15D disk revealed by 
this investigation. 



2. Observations 



2.1. Spectral Data 



Over the past few years, a number of high-resolution echelle spectra of KH 15D have 
been obtained with Star A at different elevations with respect to the grain disk edge. Some 
of these spectra include the Na I D doublet at 5895.9242 and 5889.9510 A, which allows us 
to search for evidence of absorption by neutral sodium gas which may be associated with 
the disk. These spectra were taken with the High Resolution Echelle Spectrometer (HIRES; 
Vogt et al.lll994j ) on Keck I at the W. M. Keck Observat ory in Hawaii, the Ultraviolet and 
Visual Echelle Spectrograph (UVES; iDekker et al.l 120001 ) on the Kueyen Telescope at the 
European Southern Observatory's Very Large Telescope (VLT) fac ility on Cerro Paranal 
in Chile, and the High Resolution S pectrograph (HRS; iTullI Il998l ) on the Hobby- Eberly 
Telescope (HET; iRamsey et al.l Il998l ) at McDonald Observatory in Texas. A log of the 
spectra analyzed in this paper is given in Table [TJ 

The Keck/HIRES data were taken at a resolution R = A/AA ~ 45, 000 in 2003 February 
and 2005 February (W. Herbst, P.I.) The reader is referred to lHamilton et al.l ( lin prep.l ) 
for a full desciption of the reduction procedures employed. The VLT/UVES data were 
taken at a resolution R ~ 44, 000 in two separate observing runs, one in 2001 November 
(C. Bailer- Jones, P.I.) and one in 20 04 December (R . Mundt, P.I .) Data were reduced usin g 



the UVES pipeline as described in lHamilton et al.l (120031 ) and lHamilton et al.l ( lin prep. 
The HET/HRS data were taken at a resolution R ~ 30,000 in 2005 December and 2006 



February (W. Herbst, P.I.) Data were reduced using IRAF (jTodylll993l ) routines to subtract 
the bias, perform flat-field corrections, remove scattered light and cosmic rays, extract the 
echelle science and sky orders, and calibrate the wavelength solution using images of a Th- 
Ar hollow cathode taken the same night as the science data. IDL routines were written to 
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subtract the sky spectra, normalize the science spectra, and convert to heliocentric velocities. 



Th e comparison spectr um used here is of HD 36003, which has a spectral type of 
K5 V ( ICenarro et al.l 120071 ). The spectrum was taken on the night of 2003 February 8 
with Keck/HIRES. This star has a parallax of 0.077" (jPerryman et al.lll997l ). making it 
close enough to alleviate worries about significant ISM contamination. While the spectral 
type is slightly earlier than that of the visible component of KH 15D (K6/K7), the match 
with other photospheric lines and with the extreme wings of the Na I D lines is excellent. 
The slightly later spectral class of KH 15D is offset by its lower surface gravity relative to a 
main sequence K6/K7 star. 

We initially attempted to use other comparison stars that more closely match KH 15D's 
spectral type and evolutionary phase, but HD 36003 had the best overall fit to the spectra 
of KH 15D. Spectra of Hubble 4 and LkCa 7, which are nominally T Tauri stars (TTS) of 
spectral class K7, were compared with KH 15D, but the fit with the clearly uncontaminated 
spectral features was significantly worse for both T Tauri spectra than for HD 36003. This 
is most likely due to the larger rotational broadening of the fe atures in these two T TS as 
compared to KH 15D, which has v sin i = 6.9 ± 0.3 km/s ( [Hamilton et al.ll2005l ). and 
HD 36003, whic h has v sin % < 3 km/ s ( ILuck fe Heiterl 120061 ) . Hubble 4 has v sin % 
14.6 ± 1.7 km/s (jJohns-Krull et al. I&3), while LkCa 7 has v sin i — 11 km/s ( iGlebocki et al.l 



200ll ). In addition, LkCa 7 has been noted to be a poor comparison star due to anomalous 
colors, possibly caused by a cooler companion (iGullbring et al.lll998l ). For these reasons, we 
chose to use HD 36003 as our comparison star. 

In addition to this comparison star, the bright, nearby B2 III star HD 47887 was observed 
to place constraints on interstellar Na I absorption near KH 15D ( §3.2p . 



2.2. The Na I D Profile of KH 15D 

Figure [TJ shows each of the analyzed KH 15D spectra in the region of the Na I D lines 
overlaid on the comparison spectrum, which has been appropriately shifted for each night of 
data. To accurately align the comparison star spectrum in radial velocity with the KH 15D 
spectra, we used a x 2 minimization technique based on the portion of the spectrum longward 
of the Na I Dx line (~5900-5930 A), which contains several stellar absorption features due 
to Fe I and Ti I. Formally, the comparison star spectrum should be shifted by steps in 
log(A) iTonrv fe Davisl 19791 ). but beause our analysis region includes such a small range of 



wavelengths (~15 A), the x 2 analysis is performed using uniform steps in A. The part of the 
spectrum within the wings of the Na I D lines was excluded from the fitting, because we are 
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attempting to measure very small discrepancies in the deepest part of these lines. 

Visual inspection of the KH 15D spectra overplotted on the shifted comparison star 
spectrum clearly reveals that the two match well except at the very bottom of the Na I D 
lines. Figure |2] shows an enlarged view of the region around the Na I D lines to make 
the absorption and emission features more clearly visible. Here we find that the KH 15D 
Na I D line cores deviate from those of the comparison star spectrum, showing some potential 
filling on the blue side (though telluric emission makes this difficult to see in some panels 
of Figure [Q and excess absorption on the red side of the profile. Since this absorption 
feature appears in all 12 spectra, regardless of which telescope or reduction procedure was 
employed, we are confident that it is real. As this absorption is asymmetric and not seen 
in the comparison star spectrum, this additional sodium absorption is most likely caused 
by something other than the photosphere of Star A. We turn now to an examination of the 
source of this absorption and, in particular, whether it could be associated with the gaseous 
component of the circumbinary disk. 



3. Analysis 
3.1. Measuring the Na I D Profile 



Following the procedure of iRedfieldl (120071 ). we began the process of quantifying the 



excess Na I absorption in KH 15D by taking the ratio between each available spectrum 
and the shifted comparison star spectrum after the comparison star spectrum had been 
resampled to have the same pixel sizes as the KH 15D spectrum by rebinning. The VLT 
and HET spectra are both at lower resolutions than the Keck spectra, so the comparison 
spectrum was only resampled to larger pixel sizes. The results of this process are shown in 
Figures [3] and H] where we have effectively taken out the baseline stellar spectrum, leaving 
only the excess absorption and emission. We discuss each of these in turn in § §3.1.11 and 

El 



3.1.1. Na I Emission 

While the dominant feature in these ratioed spectra is the red-shifted excess absorption 
feature clearly visible in the original spectra, this more detailed look also reveals a more 
subtle feature. All of the spectra show evidence of less absorption (or excess emission) at 
or close to the radial velocity of Star A. This emission line is especially visible in the 2003 
February Keck/HIRES data (see Figure [3]) where the radial velocity of Star A is closest to 
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the systemic velocity of 18.6 ± 1.3 km/s (jWinn et al.l 120061 ). placing the emission feature 



close to the deepest part of the absorption feature in the ratioed spectra. Regardless of the 
source of this emission, in order to accurately measure any excess absorption, this emission 
feature must be taken into account in our model. 

It is possible that the apparent excess emission is due to the spectral mismatch between 
the comparison star and KH 15D, but we believe that it is likely that Star A actually does 
have an emission core within the stellar Na I D absorption features. Three arguments favor 
the view that the excess emission is real and associated with the star. We discuss each of 
these in turn. 

First, we note that the emission feature is rather variable. On some nights the emission 
can be seen very strongly (e.g. 2005 December 22) and on others it is barely present (e.g. 
2003 February 8). While we acknowledge the danger of over-interpreting ratioed spectra, 
especially deep in the cores of strong lines, the degree of variation does seem larger than can 
easily be accounted for by simply random noise. If the apparent emission arose only from a 
mismatch of the spectral types then it should be constant in time. 

Another argument is that the emission closely matches the predicted radial velocity of 
Star A in each spectrum, even as Star A's radial velocity varies from -0.27-11.17 km/s during 
the course of our observations. This is just what is expected for emission by gas accreting 
onto the star. The predicted radial velocities of Star A are marked in each frame of Figure [3] 
for reference. 



Finally, Star A is known to have an active magnetosphere and bipolar jet ( iHamilton et al. 



20031 ). In this sense it is like a Classical T Tauri Star (CTTS), although the mass accretion 



rate is probably about an order of magnitude less than for a typical CTTS (IMundt et al. 



2010|). Emission in Na I related to the active magnetospheres and jets has been observed 



in other T Tauri stars, particularly when they are oriented such that we view them nearly 



edge-on to their disks (lAppenzeller et al.l 120051 ). as is the case here. Furthermore, in some 



spectra the em ission lines appear exten ded to the blue. This is also a common feature of 



CTTS spectra (lAppenzeller et al.l 120051 ) and it may be that we are detecting a rather weak 



version of the same phenomenon here. 

On the whole, it appears reasonable to suppose that the emission lines seen at the 
radial velocity of Star A in all of the spectra in Figure [3] are, in fact, arising within the 
magnetosphere of Star A, although the physical nature of the excess emission is not critical 
to the process of correcting for it in the calculation of the absorption column density. 

We note that in some spectra (e.g. the Keck spectra from 2003 February) there is also 
a considerable amount of night sky emission which was not fully accounted for during the 
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reduction process. Fortunately, in each case it is sufficiently shifted in radial velocity with 
respect to KH 15D that we have not needed to remove it from the line profiles because it 
does not affect the fitting results. 



3.1.2. Quantifying the Excess Na I Absorption 



In order to most accurately quantify the redshifted excess absorption, we have found it 
necessary to take into account the variable emission that follows the stellar velocity. Since the 
Na I D line is a doublet we have two absorption and two emission features to model in each 
spectrum. We have experimented with a variety of different fitting techniques employing 
constrained and unconstrained Gaussian fits to the emission and absorption components. 
Given the noise level of the data and uncertainties about the physical sourc e of the feature s 
we have, in the end, adopted a procedure quite similar to that employed by iRedfieldl ( 120071 ). 
Its implementation here is described below. 

For the two Na I D lines in the ratioed spectra, we fit a single Gaussian absorption 
component and a single Gaussian emission component to each. These four Gaussians (two 
emission and two absorption) are fit simultaneously, after being convolved with the appro- 
priate line spread function for each instrument. 

The two emission features are fixed at the predicted radial velocity of Star A for each 
night of data, with the strength varying freely. However, the area contained within the 
Gaussian curve used to model the Na I D2 emission line is fixed to be twice that of the 
Gaussian for t he Na I D\ line, due to the relative oscillator strengths of the two transitions 
(IMortonl 120031 ) . assuming optically thin Na I gas. 



The three important parameters for fitting the Gaussian absorption features are the 
column density, the radial velocity of the maximum absorption, and the Doppler width. 
Doppler width, also called the b-parameter or velocity-spread parameter, is the dispersion 
in velocities multiplied by a/2. This measure of the width of the lin e is believed to be due 
to the Doppler effect resulting from turbulent and thermal motions (jSpitzerl 119781 ) . 



For the absorption fits, the equivalent widths of the two Gaussian features were also 
fixed to be 2:1, again corresponding to an assumption of optical thinness. For our initial fit, 
the radial velocities and widths of the two Gaussian absorption features were allowed to vary 
freely and were chosen by y 2 minimization, with the Doppler widths of the two absorption 
features fixed to be the same, and the radial velocities fixed to be identical. This yielded 
absorption features with an average r adial velocity of 16 .2 ± 3.3 km/s, consistent with the 
systemic velocity of 18.6 ± 1.3 km/s (IWinn et al.ll2006l ). and an average Doppler width of 
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5.6 ± 2.5 km/s. 

For our final absorption measurements, the two Gaussian emission features at the radial 
velocity of Star A and the two Gaussian absorption features were simultaneously fit to the 
observed profiles with the Doppler widths fixed to a width that produced the best fit by eye, 
only varying the column densities (fixed 2:1) and radial velocities (fixed to be identical) of 
the lines for the x 2 minimization. This results in a measurement of the absorbing column 
density of Na I necessary to create the observed features. The results of our fits are listed in 
Table [TJ and shown in Figure EJ 

For internal consistency in our analysis when comparing spectra from different telescopes 
and instruments, error bars on the individual spectral datapoints were calculated using the 
standard deviation of a range of datapoints within the normalized continuum, so that every 
datapoint within a given spectrum is given the same error value. The error bars on the 
column densities result from the standard deviation between three methods of fitting: 1) 
freely varying all three parameters (radial velocity, Doppler width, and column density), 2) 
fixing the Doppler widths, and 3) fixing the radial velocity to match the systemic velocity. 
This is obviously only a crude estimate of the error bars, since they are estimated from 
only three measured values. However, propagation of errors results in rather similar error 
bars and the agreement among the various measurements on different nights gives added 
confidence that these error bars are reasonable. 

One source of error here is the subtraction of telluric Na I D emission lines. In the spectra 
where the radial velocities of telluric emission lines are closer to the ~18.6 km/s systemic 
velocity (e.g. the 2001 November spectrum), it is more difficult to measure the column 
densities. Another source of error is the widely varying strength of the stellar emission lines. 

An additional source of systematic error may be our slightly imperfect comparison star. 
Due perhaps to the evolutionary state of KH 15D, as discussed in §3.11 an earlier spectral 
type main sequence star actually matches the spectrum better than an identical spectral type 
T Tauri star spectrum, particularly in the deep Na I D lines. However, as noted above, using 
a non-identical comparison star could introduce a small systematic error in the derived Na I 



colum n densities, as the Na I D lines are slightly different in the K5 V star ( ITripicchio et al. 



19971 ). Such a mismatch is a possible explanation for the apparent excess flux shortward 
of Star A's central velocity, mentioned above as also perhaps arising in the magnetospheric 
emission of Star A. To a certain extent, our fitting procedure should remove this effect no 
matter what its cause, but there may well be a residual small effect on the column densities 
that is systematic in nature. 
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3.2. Correcting for Interstellar Absorption 

Figures CD and [2] show that there is indeed excess Na I absorption present for both of 
the Na I D lines. Additionally, the other absorption lines in the stellar spectrum match up 
very well with absorption lines in the spectrum of the comparison star. This alone is strong 
evidence that there is some absorbing sodium gas toward the KH 15D system. Whether or 
not this absorbing material is circumstellar needs further investigation, and in particular, 
consideration of whether the absorption could arise in the general interstellar medium (ISM) 
between us and the star. 

The only measurement of interstellar sodium in the general direction of KH 15D that 
we could find in the literature was for S Mon, a well known spectral type 07 member of 
NGC 2264 l ocated appro x imate ly 25' from KH 15D and likely at about the same distance 



from Earth. I Welty et al.l (Il994i ) measure an absorbing component at roughly the systemic 
velocity of KH 15D, but much weaker (log N^ai ~ 11.8 cm -2 ), with even weaker absorption 
extending blueward to a radial velocity of about km/s. Not only is this feature not 
strong enough to account for our observations, but the general shape does not match our 
observations. However, the distance between the sight line to this star and to KH 15D could 
result in significantly different interstellar absorption profiles. 

To get a better constraint on any possible interstellar medium contribution to the Na I D 
feature in KH 15D, we observed the bright B2 III star HD 47887, a cluster member located 
only 37" away. A spectrum was obtained with Keck/HIRES on 2003 February 9 to search 
for evidence of Na I absorption in that star. Figure [6] shows the Na I D profile of the star 
and, indeed, two components of interstellar Na I are clearly visible. We have fit them using 
the same methods described in §3] above for KH 15D, shown in the figure. The narrower, 
deeper component has a column density of log iVNai = 12.002 ± 0.002 cm -2 , with a Doppler 
width of 1.35 ± 0.01 km/s and a radial velocity of 6.24 ± 0.05 km/s according to this fit. 
The wider, shallower feature has a log N Nal = 12.093 ± 0.001 cm~ 2 , with a Doppler width 
of 10.88 ± 0.02 km/s and a radial velocity of 15.65 ± 0.05 km/s. 

This analysis shows that there may well be a significant contribution to the excess 
absorption observed in KH 15D from the general interstellar medium. However, it is unlikely 
that it can fully explain the excess absorption for the following reasons. First, the column 
densities of Na I measured in the KH 15D spectra are, in most cases, significantly larger than 
the total column density measured for HD 47887. Second, only the broad shallow feature is 
likely to be contaminating our measurements because the narrow deep feature lies at a radial 
velocity close to that of Star A in most of our spectra. Its influence would be to reduce the 
apparent amount of emission in the Na I emission feature discussed in §3.1.11 but it cannot 
be masquerading as the obvious redshifted absorption component seen in Figure [T] because 
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of the radial velocity difference. 

The broad, shallow feature, on the other hand, does nearly match the measured radial 
velocities of the excess absorption in the KH 15D spectra, but is not as deep. It is probably 
a blend of several weak components at slightly different radial velocities leading to its broad, 
shallow nature. If this absorption feature appears at exactly the same strength in the direc- 
tion of KH 15D as it does in HD 47887 only 37" away, then it is not strong enough to account 
for all of the absorption, although it does need to be corrected for in assessing how much Na I 
gas can reasonably be associated with the KH 15D disk. Our measurements have an average 
log A^Nai of 12.5 cm -2 . Subtracting the contribution of the ISM as measured in the spectrum 
of HD 47887 leaves log A^Nai ~ 12.3 cm -2 , which is our best estimate of the column density 
associated with the disk of KH 15D. Formally, the error bar is log iVNai ~ 0.1 cm" 2 , but we 
note that the possibility of systematic errors and of contributions to the column from gas 
local to NGC 2264 but not associated with the KH 15D disk means that our measurement 
may be conservatively viewed as a rather strict upper limit. Certainly the column density 
of Na I associated with the putative gaseous disk of KH 15D cannot have a column density 
much in excess of our measurement of log A^Nai ~ 12.3 cm -2 . 



4. Scale Height of the Gaseous Disk 
4.1. Comparing our Measurements to Theory 



Figure |5] shows that we have not been able to detect a significant variation of Na I 
column density as Star A moved relative to the grain disk edge. This is unfortunate because 
such a detection would have been clear evidence that the excess Na I absorption is arising 
within the gaseous component of the disk. We now inquire whether the failure to detect 
such variation is an argument that the excess absorption that we have detected could not, 
in fact, be associated with the disk. 

The exact geometry of the KH 15D system is not known and, in particular, we do 
not know the location of the sharp edge caused by the particulate disk. Our line of sight 
may intercept the disk (it is perhaps more accurately referred to as a ring in this context) 
at its inner edge, outer edge, or so mewhere between (these possible disk geometries are 
illustrated very clearly in Figure 6 of IWinn et al. 120061). The inner edge is possibly at the 
3:1 resonance with the binary system (~0.6 AU; iHerbst et al.l 120081) and the o uter edge is 
probably limited to 5 AU or less in order to precess ( jchiang &: Murray- Clay 2004 ). Assuming 
that the particulate disk marks the mid-plane of the gaseous disk and that our line of sight 
is nearly along the disk plane, we may then equate the apparent distance of Star A from 
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the opaque disk edge (A) with the height of the gas layer that our line of sight to Star A 
is probing. Of course, this is a simplification because the disk must be warped and because 
we are not really looking exactly along it. In the absence of more detailed knowledge of the 
system's geometry, it is a reasonable approximation. 

The scale height (H) of the gas layer may be estimated under the assumption of hydro- 
static equilibrium in the direction perpendicular to the disk (z). Assuming an ideal gas at a 
temperature (T) and constant vertical component of the local gravity (g z ), one may write 



H = g 



-i 



g; 1 



kT 
LfirriH 



where k is the Boltzmann constant, /i is the mean molecular weight of the gas and m# is the 
mass of the hydrogen atom. We further assume that the disk gas has a cosmic abundance of 
helium, is composed of predominantly neutral atoms, and that the Na atoms are well mixed 
with the H and He. In this case, we may set fi = 1.3. 

T is expected to be in the range of about 150 to 400 K, depending on the exact location 
of the edge. The local vertical component of gravity, g z , may be set by either the self-gravity 
of the disk or by the vertical component of the stellar gravity. In the first case, g z = 27rGS, 
where £ is the mass per unit area in the local disk, and in the second case q 7 = Q 2 z for 



small z, where f2 is the angular velocity of the disk material (IGoldreich fc Ward! 1 19731 ). For 
us, the more interesting limit is the weakest gravitational field that could be present, so 
we ignore any possible self-gravity. The range of motion of Star A in our data set is only 
about 3 stellar radii (Ra, where Ra = 1-3 -R©), so we take z = 1 (in units of Ra) for this 
approximate calculation. 

With the assumptions above, we estimate that a putative gas disk in hydrostatic equi- 
librium would have a measured scale height of H = 14 (2600) Ra at the inner (outer) edge 
of the disk, located at 0.6 (5) AU. These scale heights are shown in Figure [5j This is too 
large a scale height for us to have measured, given the limited range in A probed by Star A 
during our observations and the precision of our measurements. Conversely, one can argue 
that if the the slightly decreasing trend of Na I column density with A, visible in Figure |5l 
is real, then the occulting edge of the gas disk is likely to be the inner edge of the disk and 
there may well be a contribution to g 2 from the self-gravity of the disk (i.e. it would need to 
be a few times larger than the stars alone could provide). 
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4.2. The Settled Grain Disk of KH 15D 



We did not detect any statistically significant dependence of the excess Na I absorption 
on the height of Star A above the occulting edge. One expects from the order-of-magnitude 
analysis in §4. II that the scale height will be much greater than the few stellar radii range of 
A for which we have data. By contrast, the scale height of the grain disk responsible for the 
extinction must be less t han 1 Ra- It has been represented in m odels based on the light curve 



as a sharp "knife edge" (iHerbst et al.l 



2002 



Winn et al.l 120061 ) . and may be estimated from 



Figure 1 of the Supplement of lHerbst et al.l ( 120081 ) as having the value 0.7 Ra- While there is 
uncertainty in all of these estimates it seems clear that the scale height of the particulate disk 
is at least one order of magnitude less than the scale height of any gaseous disk associated 
with it. 

Furthermore, our lack of detection of a strong Na I feature associated with the steep 
increase in extinction at the grain disk boundary demonstrates again how different (evolved) 
this obscuring matter is from interstellar grains. The extinction in the I band for our most 
obscured data point (the first one in Table [1]) is 0.9 mag, which for a normal interstellar 
medium extinction law corresponds to 1.6 mag of extinction in the V band. The reddening 
associated with this amount of extinction would be E(B-V) ~ 0.5 m ag, which would predict 
a Na I column density of log A^f a i — 13.2 cm -2 (IBohlin et al.lll978l ). much higher than the 
Na I column density that we actually detect. 

Our interpretation of these facts is that the particulate matter in the KH 15D disk has 
grown large enough to become both optically grey and dynamically separated from the gas 
layer. This req uires particle sizes o f a millimeter or so, consistent with the minimum grain 
size derived by lHerbst et al.l (120081 ) on the basis of a completely different argument. 



It would be interesting to know what the gas-to-dust ratio in the settled particulate 
layer is. Unfortunately, converting from Na I column density to total gas column density 
is difficult for two reasons. First, the ionization fraction of sodium atoms is unknown, and 
second, it is likely that much of the sodium is depleted onto the solids. If ISM values of 
the ionization and depletion applied, then there would be significantly more matter in the 
form of solids than gas. If the gas density is very high, on the other hand, one would expect 
the solid grains to quickly lose momentum through friction with the gas and spiral into the 
stars. Unless there were a steady source of new dust, such a structure would disappear 
on an astronomically brief time scale. It may be, therefore, that the very existence of the 
particulate disk indicates a fairly low gas to dust ratio. 
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5. Summary and Conclusions 



Na I lines in the spectrum of the binary WTTS KH 15D have been analyzed in detail. 
We find an excess absorption component that may be attributed to the gas component of a 
circumbinary disk. The derived column density is log iV NaI ~ 12.5 cm -2 with no significant 
variation associated with the position of Star A relative to the optically occulting edge. 
Subtracting the likely contribution of the ISM leaves log iVNai ~ 12.3 cm" 2 . 

There is no detectable change in the gas column density across the "knife edge" formed 
by the edge of the grain disk, indicating that the gas and solids have very different scale 
heights, with the solids being highly settled. 

If a standard ISM ratio of A^ai to Ay is applied along these lines of sight, there would 
be a much lower mass of gas than there is mass of solids in this disk. However, the conversion 
from Na I column density to total gas column density is complicated by the fact that much 
of the Na in this disk gas might be either ionized or depleted onto the solids. 

An additional complication is detectable emission in the Na I feature that follows the 
radial velocity of Star A. While this might be due to the slight mismatch in the spectral type 
of the comparison spectrum, it more likely arises from the active, accreting magnetosphere 
known to be associated with this (quasi) WTTS. 

Our data support a picture of this circumbinary disk as being composed of a very thin 
particulate grain layer composed of sand (millimeter) sized objects or larger, which have 
settled within whatever remaining gas may be present. This phase of disk evolution has 
been hypothesized to exist as a prelude to the formation of planetesimal s through gravita- 
tional fragmentation (jSafronov &: Zvjaginalll969l ; iGoldreich fc Wardlll973l ) but has not been 
extensively observed, and is expected to be short-lived if much gas were still present in such 
a disk. 



Both compo nents of KH 15D are now fu lly eclipsed by the circumbinary disk for the 
forseeable future (jChiang fc Murray-Clayll2004l ). so our chance to directly measure absorption 
in this particul ar disk has passed. F ortunately, other young, eclipsing, edge-on disk systems 
such as WL4 (IPlavchan et al.l 120081 ) and YLW 16A (Plavchan et al., in prep.) are being 
discovered in photometric monitoring surveys, and future spectroscopic studies may allow 
better understanding of the timescale of dust and gas settling in transition disks. 
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Table 1. Na I Column Densities 



Observation 
Julian Date 


Observation 
UT Date 


Observatory/ 
Instrument 


(Ra) 


Star A RV 6 
(km/s) 


BC C 
(km/s) 


log iV Na i 
(cm -2 ) 


2453354.7199 


2004 Dec 15 


VLT/UVES 


-0.43 


7.53 


8.73 


12 


60 


± 





17 


2452679.8276 


2003 Feb 9 


Keck/HIRES 


-0.39 


11.17 


-18.71 


12 


63 


± 





07 


2453429.8250 


2005 Feb 28 


Keck/HIRES 


-0.07 


4.09 


-25.18 


12 


31 


± 





18 


2453353.6774 


2004 Dec 14 


VLT/UVES 


0.52 


6.29 


9.29 


12 


75 


± 





20 


2452678.8193 


2003 Feb 8 


Keck/HIRES 


0.75 


9.53 


-18.29 


12 


40 


± 





03 


2453430.8201 


2005 Mar 1 


Keck/HIRES 


0.92 


2.89 


-25.45 


12 


38 


± 





15 


2453352.7976 


2004 Dec 13 


VLT/UVES 


1.25 


5.35 


9.56 


12 


82 


± 





48 


2453724.7264 


2005 Dec 20 


HET/HRS 


1.57 


0.33 


6.35 


12 


05 


± 





32 


2453726.7222 


2005 Dec 22 


HET/HRS 


2.02 


-0.18 


5.35 


12 


73 


± 





34 


2453780.7305 


2006 Feb 14 


HET/HRS 


2.11 


0.38 


-20.45 


12 


66 


± 





96 


2453727.7181 


2005 Dec 23 


HET/HRS 


2.16 


-0.27 


4.84 


12 


45 


± 





07 


2452242.7446 


2001 Nov 29 


VLT/UVES 


3.12 


8.60 


16.23 


12 


08 


± 


1 


04 



a A is the height of the center of Star A above the grain disk edge, measured in units of 
the radius of Star A {Ra = 1.3 R & = 9.0 x 10 10 cm). Note that because A is measured in 
stellar radii, none of these spectra were taken when Star A was completely obscured by the 
disk. 



b Calculated based on a slightly refined version of Model 3 of IWinn et all ( 120061 ). 



c Barycenter correction, which shows the radial velocity telluric Na I features will have in 
each of our spectra. 
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Fig. 1. — Spectra of KH 15D from different nights and the comparison star spectrum, which 
has been shifted to match the KH 15D spectra. The KH 15D spectra are shown in black, 
while the comparison star spectrum is red. Nights are arranged in order of increasing A (see 
Tabled]). Note that absorption features are very well aligned, and that the Na I D lines both 
show obvious excess absorption on the positive radial velocity side. 
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Fig. 2. — Zoomed in view of the spectrum from Keck taken on 2003 February 8 (black) 
and the properly shifted comparison star spectrum (red), allowing closer inspection of the 
absorption and emission features. 
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Fig. 3. — KH 15D spectra from each night ratioed with the comparison star spectrum (black), 
shown in order of increasing A (see Tabled]). Also shown are the Gaussian absorption feature 
fits to the excess Na I absorption (blue filled) and the modeled emission plus absorption (red). 
Vertical lines mark the radial velocities of telluric Na I emission (dotted), Star A (dashed), 
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Fig. 4. — Zoomed in view of the spectrum from Keck taken on 2003 February 8 ratioed with 
the comparison star spectrum (black), allowing closer inspection of the Gaussian absorption 
feature (blue filled) and the modeled emission plus absorption (red). Vertical lines follow 
Figure [3l 
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Fig. 5. — The measured column densities plotted against A, the elevation of Star A above 
the grain disk (Keck: blue diamonds, VLT: green triangles, HET: red squares). Only two 
HET/HRS datapoints are shown, because the three datapoints with A close to 2 Ra were 
averaged together, weighted by their error bars, to one point. The error bars show the 
systematic errors, as discussed in §3.1.21 The three orange lines show 3 different scale 
heights: the scale height of the grain disk (H = 0.7 Ra', dotted), the scale height calculated 
at the inner edge of the gas disk (H = 14 Ra; solid), the scale height calculated at the outer 
edge of the gas disk (H = 2600 Ra', dashed). See §4.11 for calculations. 
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Fig. 6. — Interstellar Na I absorption toward the spectral type B2III star HD 47887. Two 
absorbing components are seen. The blue filled areas show the separate Gaussian components 
before convolution with the line spread function, while the red solid line shows the sum of 
the two components after convolution. The solid vertical line shows the systemic velocity of 
KH 15D for reference. 



